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In a recent letter [l[ , we demonstrated that there exists 
a nontrivial correction, arising from the viscosity of the 
electron liquid, to the conductance of nanoscale junctions 
calculated within the adiabatic local-density approxima- 
tion (ALDA) to time-dependent density-functional the- 
ory (DFT). This dynamical correction cannot be cap- 
tured by any static DFT functional, even the exact one. 
We showed that the effect of the viscosity on the conduc- 
tance can be analytically predicted, in a qualitative way, 
using time-dependent current-density functional theory 
in the zero-frequency and linear-response regime [3] ■ In- 
deed, in the DC hmit, we found that these viscous effects 
increase the resistance. In order to provide an estimate 
of these effects, we derived Eq. (14) for the dynamical 
viscous resistance, and evaluated that expression by em- 
ploying the viscosity of bulk electrodes. 

Jung et al. have used our Eq. (14) to calculate the dy- 
namical resistance of two infinite jellium electrodes sepa- 
rated by a vacuum gap ^ . By using the local-density de- 
pendent formula for the viscosity coefficient as reported 
in Ref. Q they conclude that the dynamical resistance is 
negligibly small in the specific cases they consider. Here 
we argue that their calculations do not preclude the pos- 
sibility that the viscosity contribution to the conductance 
be large in realistic nanoscale structures, the systems of 
interest in Ref. fH. Two main points support this state- 
ment. 

1) Our Eq. (14) is an approximate formula, derived 
under certain physical assumptions (see below) in order 
to provide a qualitative understanding of the viscous ef- 
fects in nanostructures. As noted in Ref. [l|, Eq. (14) 
was derived assuming homogeneous density in the trans- 
verse direction and homogeneous current density in both 
transverse and longitudinal directions. The only contri- 
bution to the correction we included in Eq. (14) comes 
from density variations along the longitudinal direction. 
It is on the ground of these assumptions that we adopted 
the viscosity of the bulk electrodes in the model calcula- 
tions. We did not claim any quantitative accuracy of our 
estimates in Ref. [l[. 

In realistic nanoscale structures, the current density 
may vary rapidly in the transverse direction due to a de- 
crease of velocity from the center of the channel to the 
sides of the conductors. The transverse density and cur- 
rent density gradients can thus contribute significantly 
to the dissipation effects. The contribution can be fur- 
ther enhanced in the presence of turbulent eddies near 
the contacts [1, To capture these gradient contri- 
butions, one needs to evaluate the dissipation power 
dE/dt — — J ej ■ Excdr directly (and the associated resis- 
tance as R'^y'' = {dE/dt)/ P, with I the total current), 
because the nonconservative nature of the dynamical xc 
field makes it, in general, ambiguous to evaluate a line 
integral as in Ref. [l|. If, once again, the current density 
and viscosity are assumed constant, the correction to the 
resistance evaluated from the power dissipated is given by 



Rdyr^ = / [4(5^„-i)2 + (a^„-i)2] dr, where ± rep- 
resents the transverse direction, rj the (constant) viscos- 
ity, and A is the cross section of the nanostructure. This 
expression contains a positive transverse density gradi- 
ent term that has been neglected in Eq. (14) of Ref. [l[. 
This term thus increases the dynamical effects evaluated 
in Ref. [ij. More generally, the transverse density and 
current density gradients and the spatial variation of the 
viscosity must all be taken into account when evaluat- 
ing the viscous resistance. For a general current-carrying 
nanoscale system, therefore, a quantitative evaluation of 
these corrections requires knowledge of their microscopic 
current and density distributions [5j , and the dissipation 
power (and associated resistance) must be evaluated nu- 
merically. 

2) To calculate the viscosity coefficient, Jung et al. 
have applied the Conti-Vignale formula, Eq. (4.10) of 
Ref. This formula is more accurate than the high- 
density formula we have used in [l|. However, unlike 
what Jung el al. suggest, the Conti-Vignale formula does 
not interpolate between the high-density and low-density 
limits of the homogeneous electron liquid. The simple 
reason for this is that the exact low-density limit of the 
viscosity of the electron gas is unknown. The formula in- 
stead comes from a fit to numerical results of Nifosf et al., 
which are based on mode- mode coupling theory Q • This 
theory is certainly not exact in the low-density limit 
Indeed, due to the strong correlation effects in the low- 
density electron gas, where the electrons are on the verge 
of crystallization, it is reasonable to suspect that the rel- 
ative viscosity rj/n might increase well above Eq. (4.10) 
of [3| with decreasing density. 

We conclude that the comment of Jung et al., while 
interesting, should not be taken as an indication that the 
viscosity corrections to the conductance of real nanoscale 
structures are small. A more accurate treatment of the 
density and current density distribution and of the elec- 
tronic correlations may yield much larger corrections in 
realistic systems. 
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